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Abstract 


We  discuss  the  decay  of  a  system  of  currents  in  a  (bounded)  region 
of  finite  conductivity,  assuming  that  the  decay  is  very  slow. 

The  chief  results  are  a  set  of  inter-related  (vector)  extremum 
principles  for  calculating  the  decay  times  associated  with  regions  of  arb- 
tray  shape  and  conductivity.  These  expressions  differ  according  as  the 
electric  field  vector,  the  magnetic  field  vector,  or  both,  are  permitted 
arbitrary  variations. 

The  decay  times  are  shown  to  be  bounded  from  above;  therefore,  the 
values  obtained  from  the  variational  principles  aoproach  the  exact  decay 
times  from  below,  and  the  stationary  property  corresponds  to  the  theorem 
that  the  exact  field  makes  the  decay  time  a  maximum  (Lena's  law), 

A  vector  iterational  procedure,  generalizing  the  standard  procedure 
in  scalar  problems,  is  discussed,  and  it  is  pointed  out  that  two-dimensional 
decay  problems  also  may  be  treated  by  conformal  mapping, 

A  numerical  application  is  the  subject  of  a  companion  report* 
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1,  Introduction 

Suppose  a  system  of  ctorrents  has  been  established  in  a  certain 
region.  If  the  conductivity  is  everywhere  either  zero  or  infinite,  the 
fields  persist  indefinitely.  However,  If  regions  of  finite  conductivity 
are  present,  they  eventually  convert  all  the  field  energy  into  joulean 
heat. 

We  wish,  in  what  follows,  to  consider  such  magnetic  decay  processes 
in  the  limit  of  large  decay  times.  To  make  this  restriction  more  precise, 
we  first  note  that,  if  c  is  the  velocity  of  light,  ^  the  decay  time,  a    the 
conductivity  (e.s.u. ),  and  L  a  typical  dimension  of  the  region,  then,  as 
shown  later, 

(1.1)  r-(y-4    • 

c 

Using  (1»1),  we  may  state  our  limiting  condition  in  two  equivalent  forms: 
We  shall  restrict  our  discussion  to  situations  for  which 

(1.2)  (cT/L)^  '^^  crT»  1   • 

Physically,  the  first  expression  states  that  the  decay  times  are  very  large 
compared  with  the  time  an  electromeignetic  disturbance  requires  to  traverse 
the  region^  Lhe  second  expression,  applicable  only  where  the  conductivity  is 
different  from  zero,  states  that  the  conduction  current  is  much  larger  than 
the  displacement  current.  Mathematically,  (1.2)  means  that  the  simplified 
equations  of  the  quasi- steady  theory  (no  displacement  cuirrent)  are  applicable. 
Our  interest  in  this  problem  stems  from  recent  estimates,  by  Cowling'-  ■' 
and  Wrubel"-  -^ ,  of  the  decay  times  associated  with  stellar  magnetic  fields. 
Cowling,  using  a  variational  method,  calculated  the  decay  time  characteristic 
of  the  lowest  mode  of  the  solar  magnetic  field,   Wrubel,  using  n\imerical  in- 
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tegration,  confuted  decay  times  for  several  higher  modes.  The  present  remarks 
start  from  the  observation  that  Cowling's  variational  method  may  be  generalized 
to  apply  to  an  arbitrary  vector  mode,  or  indeed  to  an  arbitrary  vector  field. 

In  what  follows,  we  give  a  systematic  development  of  a  vector  formalism 
suitable  for  the  calculation  of  decay  times  associated  with  regions  of  arbitrary 
shape  and  conductivity.  Our  chief  results  are  a  set  of  interrelated  extremum 
principles  for  the  decay  time.  These  expressions  differ  from  one  another  accord- 
ing as  the  electric  field  vector,  the  magnetic  field  vector,  or  both,  are  per- 
mitted to  undergo  arbitrary  variations  about  their  correct  values.  Since,  as 
we  shall  show,  the  decay  times  are  bounded  from  above,  our  variational  principles 
approach  the  exact  decay  times  from  below.  This  fact  may  be  expressed  as  a  theorem: 
the  exact  field  makes  the  decay  time  a  maximum. 

A  vector  iterational  procedure,  for  use  with  the  variational  principles, 
can  be  set  up  quite  simply;  such  a  scheme,  generalizing  the  standard  procedure 
in  scalar  problems,  is  described.  As  a  final  remark,  it  is  pointed  out  that  two- 
dimensional  decay  problems  may  be  treated  by  conformal  mapping  procedures. 

In  this  report,  we  give  a  general  discussion  of  the  mathematical  formalism, 
describing  the  variational  principles  and  the  iterational  procedure  in  detail.  In 
a  companion  report  •-  -' ,  which  may  be  read  independently,  one  of  the  variational  ex- 
pressions developed  here  is  derived  directly  and  is  applied  to  the  problem  consider- 
ed by  Wrubel.  Comparison  with  Wrubel's  nvonerical  integration  shows  that  a  simple 
trial  function  approximates  the  exact  result  quite  closely. 

2.  The  problem 

We  consider  a  set  of  currents  in  a  body  of  conductivity  ctCj).  The  cur- 

-t/T 
rents,  and  the  fields  associated  with  them,  decay  according  to  e  '      ,     We  wish 


to  calculate  'X » 
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We  shall  be  interested  only  in  the  quasi-steady  case.  That  Is,  (1.2) 

applies.  Then,  from  the  first  equation  below,  E  '^  H  (L/c  'J')}   in  the  second 

equation,  therefore,  the  displacement  current,  -  -—  ,  can  be  omitted,  since  it 

c~ 

is  of  order  (l/c  T)  compared  with  the  curl  H  term. 

Thus,  the  field  equations,  in  Gaussian  units,  are  (e  »  |j.  =  1): 

(2.1)  V  A  E  -  H/c  r 

(2.2)  V  X  H  -  ^  o-E 

c 

where  the  last  equation  will  also  be  valid  outside  the  body,  within  a  region 
of  dimension  •^L,  if  we  put  cr  ■  0. 
We  shall  also  require  that 

(2.3)  V  »  H  -  0 
(2.U)            V  •  E  -  0  . 

The  first  of  these  equations  is  a  consequence  of  (2.1 )j  the  second,  which 
states  that  the  charge  density  is  always  zero,  is  closely  related  to  (1.2 )i 
for,  if  o'T»  1,  the  volume  charge  decays  in  a  time  that  is  very  small 
compared  to  the  decay  time. 

A  further  condition  follows  from  Eqs.  (2.2)  and  (2.U): 

(2.5)  E  •  y  o-  -  0  . 

For  exanple,  if  (y  is  spherically  symmetric,  E  will  not  have  a  radial  component 
(H-modes  only). 

Because  Eqs.  (2.1)  and  (2.2)  are  valid  only  within  lengths  -^L,  it  is 
convenient  to  limit  our  analytical  manipulations  -do  some  region  of  maximum 
dimension  <  L.  The  conducting  body,  assumed  bounded,  itself  suggests  such  a 
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region.  The  influence  of  the  fields  outside  the  body  is  represented  by 
imposing  a  boundary  condition  on  the  surface,  S,  of  the  body.  Using  the 
(sjrrametric)  dyadic s,  ^  and  ZY,  and  the  unit  normal,  n,  pointing  out  of  the 
body,  the  boundary  condition  can  be  given  in  either  of  the  two  equivalent 
forms 


(2.6) 

on  S. 


n  X  E  -  ^  .  h1 

(2.7)  1  ^  n.  'J^    .  E 

Here,   the   (surface)  dyadics    U     and    ^    are  related  by: 

(2.8)  y  »  -n    Y  y^    )f    nj  ^  -  -n^-"-"-  x    n. 

Outside  S,  both  E  and  H  satisfy  Laplace's  equation. 

Poynting's  theorem  follows  from  Eqs.  (2.1)  and  (2.2): 

(2.9)  V  .  (£  X  H)  -  1^  -  i^  tf-f 
as  do  the  additional  identities; 

(2.10)  2l  +  iii  o"!^  »H»VxE  +  E.VX_H 

"  V»(ExH)+  2E.  •  7XH^ 
■  '•(H.''E)+  2H.  •  VxE.. 

In  the  following,  we  always  integrate  over  the  conducting  region 
(volume  V,  surface  S).  For  example,  from  (2.9), 

f  dV  H^ 

(2.11)  cr  -  1—^4 ^ J . 

i^  J  dV  0- E;^  +  J  dS  n  •  E  X  JL 

This  equation  is  a  field  generalization  of  the  definition  of  the  decay  constant 
of  a  ivunped  constant  circ\iit:  for  a  simple  R,  L  circuit,  the  decay  time  is 
'*'  L/Rj  in  oyyr   case,  (2.J.1J  defiiies  X    as  the  ratio  of  the  stored  magnetic 
energy  to  the  energy  dissipated  per  second. 
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Although  ^  may  be  computed  from  (2.1i),  the  stationary  properties 
of  this  expression  are  not  clear.  These  properties,  which  are  more  conven- 
iently based  on  the  identities  (2»10),  will  now  be  discussed, 

3»  Variational  expressions 

Eqs.  (2.1)  and  (2.2),  and  the  boundary  condition  (2.6),  may  be  derived 
by  demanding  that  the  expression  below  be  stationary  when  E  and  H  are  varied 
independently: 


(3.1) 


^  ■  I  j  H  ^  *  ^  ^^^  "  2^  •  '  ^  hI  -  I  fdS  H  .  ^  .  H 


If  the  boundary  condition  is  used  in  the  form  (2.7 )>  the  Lagrangian  above  is 
replaced  by: 


(3.2)     L^  » 


We  now  verify  the  stationary  properties  of  L.  and  Iv,; 

51^  -         dViJ  ^  cr^  *  T"  ^-  '   ^-^-^I*    VxH-E.   V     6hI    -       dS  5g  •%  •  H 
The  last  term  in  the  volume  integral  may  be  rewritten; 

(3.3)  rS  •  V  ^  5  H  -  V  .    (EX6H)  -6H  •   7  XE, 

so  that: 

(3.'-i)  61^-      j     dTi  6a  .     A-VxE     +  5E  .    l^crE-7xH| 

+      I      dS  6H  •      n>£  -^     •    H  1    . 


« 


The  Lagrangians  used  in  this  sections  were  suggested  by  analogous  forms, 
for  the  oscillatory  case,  discussed  by  J.  Schwinger  in  war-time  lectures 
at  the  M.I.T.  Radiation  Laboratory,  sind  by  N.  Marcuvitz  in  recent  lectures 
at  New  York  University. 


-  6  - 

Similarly, 


(3.5) 


Dig  -  I  dvhn  .  TA  -7  X  e|  +  6E  •  \^  o-E  'Vxm 


-   dS 


5E  •  Hxn  -U    •  e1 


We  conclude  from  these  equations  that  for  independent  variations  of 
E  and  H  about  the  values  that  satisfy  Eqs,  (2.1),  (2.2),  and  (2.6),  L.  is 
stationary^  conversely,  if  I^  is  stationary  for  independent  variations  of 
E.  and  ii,  Eqs.  (2.1),  (2.2),  and  (2.6)  follow.  A  corresponding  statement 
holds  for  L-  if  Eq.  (2.6)  is  replaced  by  (2.7).  Also,  from  the  identities 
(2.10),  the  correct  fields  -  and  boundary  conditions  -  make  the  stationary 
values  of  both  L.  and  Lp  zero.  Thus,  if  L,  and  L_  are  each  equated  to  zero, 
we  get  two  variational  expressions  for  the  decay  time,  T". 

In  these  variational  expressions,  both  E,  and  H_  can  be  varied  independently. 
As  an  alternative,  we  may  use  one  Maxwell  equation  to  define  one  of  the  field 
vectors  in  terms  of  the  other.  Putting  this  jnto  our  variational  principles 
gives  us  two  new  variational  expressions,  each  involving  only  a  single  field 
vector.  Varying  this  vector,  and  imposing  the  stationary  property,  yields  the 
other  Maxwell  equation  and  the  boundary  condition. 

For  example,  we  may  define  E  in  terms  of  H  by  Eq.  (2.2).  Putting  this 
in  L. ,  we  get  L.  t 


(3.6) 


^■^{^-^]-|/-a./.K. 


Using  (3.3),  one  easily  verifies  that  61^  »  0  for  variations  of 
H^ about  (2,1)  and  (2.6).  Since  L,  ■  0  for  the  correct  choice  of  H,  (3.6), 
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set  equal  to  zero,  is  a  variational  principle  for  ^  in  which  only  the  magnetic 
field  is  varied, 

A  similar  variational  principle  for  ^  ,  that  involves  only  the  electric 
field,  is  obtained  by  defining  H.  in  terms  of  E.  by  Eq.  (2.1);  putting  this  into 
L_j  we  find: 

h'l  J    ^^  ]  1^  <5-E^  -  cT  (V  X  E)i  +  I   f  dS  E  .  t^  •  E  . 

t  I 

It  now  follows  that  L_  »  0  if  E  satisfies  Eqs.  (2.2)  and  (2,7);  also  SL^  "  ° 

if  E  is  varied  about  the  solution  of  these  equations.  Hence,  (3.7), set  equal 
to  zero,  is  a  variational  principle  for  7"  that  involves  electric  field  varia- 
tions only. 

Using  Eqs.  (2.6)  and  (2.7),  we  note  that: 


(3.7) 


(3.8) 


n«E.AH  ■  H'J^'H  -  E  -^  •  _E. 


The  stationary  properties  of  (2.11)  are  now  clear:  E  may  be  eliminated  from 
(2.11)  by  using  (2.2)  in  the  volume  integral  and  (3.3)  in  the  surface  integral; 

the  result  is  identical  with  the  variational  principle  obtained  by  putting 

t 

L.  »  0.  Similai'ly,  H  may  be  eliininated  from  (2,11)  by  using  (2«l)  in  the 

volume  integral  and  (3»3)  in  the  surface  integral,  the  result  now  being  the 
variational  principle  obtained  by  setting  L^  ■  0. 

From  an  examination  of  (3.U)  and  (3»5)j  we  see  that  the  sole  function 
of  the  surface  terras  in  (3«1)  and  (3.2)  is  to  insure  that  the  variations  are 
made  about  fields  that  satisfy  (2.6)  and  (2.7).  Correspondingly,  in  (3.6)  the 
surface  term  limits  the  allowable  variations  of  H:  the  magnetic  field  must  be 
varied  about  the  field  that  satisfies 

(3.9)  J3.X  flj^]        .  :^  .  H  (on  S) 


w  -y^- 
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A  like  statement  holds  i'or  (3*7):  because  of  the  surface  term,  the  variations 
of  E  must  be  made  about  the  field  that  satisfies  (we  have  set^j^  »  ^'^J^   ^* 

(3.10)  (V  X  E)  X  n  'U'*   .  E  (on  S) 

It  follows,  then^that  if  we  impose  the  restrictions  (3»9)  and  (3»10)  separately 

t         I 

from  the  variational  principles,  the  surface  terms  in  L.  «  0  and  I^  •  0  may  be 

eliminated. 

For  example,  setting  L^  ■  0,  and  using  (2.2),  (3.3),  and  the  identity 


W-  '-l^h^-'Mi^S-, 


we  find: 


(3.11)  cT 


J  dV  Hf 


/<i^^-  '^^1^ 


(3.11)  is  stationary  for  variations  of  H  about  the  field  that  satisfies  {,^,,9) 
on  the  surface  of  the  conducting  body  and 

in  the  interior. 

If  we  set  L-  «  0,  we  may  \xse  (2.1)  and  (3.8)  to  get: 

fdviiirf 

(3.13)  cr  -  -^ ~ 


J  dV  E.*   VX(V  A- E) 

This  expression  is  stationary  for  variations  of  E.  about  the  field  that  satisfies 
(3.10)  on  S  and 
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(3.114)  V  A-  (V  A  E)  -  k^^ 

inside  S. 

We  conclude  this  section  with  a  brief  discussion  of  the  spectrum  of  '^« 
We  put 

(3.15)  X  =  ^ 

c^r 

and,  as  required  by  (1.2),  we  impose  (2.1^)  on  (3.lU).  The  result  is: 

(3.16)  (7^  +  X<r)  E  =  0, 

where  E  must  be  represented  by  its  rectangular  components.  We  now  wish  to 
show  that  the  spectrum  of  ^   (as  determined  by  the  differential  equation  (3.16) 
and  the  boundary  condition  (3.10)y  is  bounded  from  abovej  for,  if  this  is  so, 
our  variational  expressions  for  ^    will  always  approximate  T  from  below. 
We  shall  sketch  a  proof  based  on  quantum-mechanical  reasoning:  Put 

(3.17)  ^  '  ^<^  '':¥  V. 

Then  (3.16)  becomes,  for  each  component  of  E^,  a  Schrodinger  equation  that  des- 
cribes a  particle  of  mass  m  bound,  with  zero  binding  energy,  to  a  potential  V. 
(Since  Xo"  is  always  >  0,  V  is  an  attractive  potential.)  Thus,  the  equivalent 
quantum-mechanical  problem  concerns  a  potential  of  fixed  range  (L)  and  shape  (o")j 
the  scale  factor  of  the  potential,  that  leads  to  zero  binding  energy,  is  the 
eigenvalue  that  must  be  determined. 

Since  the  particle  energy  is  zero,  and  V  is  always  <  0,  the  magnitudes 
of  the  kinetic  energy,  T,  and  the  potential  energy,  V,  are  always  equal:  That 
is, 

IVI^T.  |!  , 
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where  P  is  the  particle  moraentvun.     On  the  other  hand,   because  of  the  uncertainty 
principle,   P  is  always  of  the  order  of,   or  greater  than,   n/L,     Therefore 

f^     |V|    l2>1 
n 

or,  returning  to  our  decay  problem  via  (3.17), 

(3.18)  ^<    ^ 

c 

■a- 
Thus,  T  is  bounded  from  above  ,  and  the  values  of  T    given  by  using  approx- 
imate fields  in  the  variational  principles  will  always  be  less  than  the  true 
value  of  T. 

These  remarks  may  be  expressed  as  a  thoerem:  the  exact  field  makes  the 
decay  time  a  maximum.  In  this  fomi,  we  have  what  appears  to  be  an  example  of 
Lenz's  law:  namely,  the  fields  adjust  themselves  to  oppose  any  change. 

Incident ly,  we  note  that  (3.18)  justifies  the  estimate  of  T   given  in 
(1.1). 

U.  Iteration  of  £ 

The  variational  expressions  given  above  are  most  useful  when  combined 
with  some  systematic  method  of  improving  the  trial  functions  to  be  used.  Such 
an  iterationai  scheme  is  well  known*-  ^   for  scalar  variational  principles  of  the 
type  we  have  discussed.  We  wish  to  show,  in  this  section,  that  an  iterationai 
procedure  may  also  be  set  up  for  the  vector  variational  principles  of  Section  3. 

A  simple  way  of  accomplishing  this  is  merely  to  generalize  the  usual 
scalar  treatment  to  apply  to  vector  functions.  Thus,  we  define  the  sequence  of 
functions  K,  by  the  conditions: 

#  —————— —^ 

The  argument  may  be  extended  to  show  that  the  spectrum  of  ^  is  also  discrete. 


-  11  - 


(U.l)  V  .  E^  -  0 

(U.2)  V  X  (V  X  Ej^)  =  -7^  ^  -  ^    <rE^^^ 

c 

ih.3)  (7  >  Ej^)  A  ji  -   Z^'  •  E^  (on  S). 


In  the  second  equation,  the  Lapiacian  operates  on  the  rectangular  components 
of  E^. 

If  we  apply  this  iterationai  scheme  to  Eq.  (3 ,13),  we  can  show  that, 
as  k  -»  00,  the  trial  functions,  Ej.  and  the  corresponding  values  of  'T^  ,  Tn^y 
approach  the  exact  eigenfunctions  and  eigenvalues* 

The  proof  follows  in  close  analogy  with  the  scalar  case:  First,  using 
(U»3),  we  note  that: 

J   dS  Ej^    *  J"        0     Ej^  -  J   dV  J  (V  A  E^)  .  (V  x  ^)  -  E^  •  V  A  (V  x  Ej^) 

The  left  side  of  this  equation,  and  the  first  term  on  the  right,  are  symmetrical 
in  Ji  and  k.  Therefore,  the  remaining  term  on  the  right  must  also  be  unaffected 
bj'-  interchanging  Z.    and  k.  That  is,  introducing  (_£  +  k  -  1)  as  defined  below. 


(U.U) 


(/  +  k  -  1)  -  I  dV  <rY^    .  E^_^  «  j  dV  o-E^_^»  \  , 


the  integral  in  (U«U)  depends  only  on  the  sum  of  the  indices. 

If  we  put  E  into  (3*13)  and  call  the  resultant  approximation  'T^,. 


we  find: 


(U.5)  r. 


(2k) 


2k    (2k-l)  * 

We  now  show  that  the  sequence  'T„,  increases  with  k.  The  proof  pro- 
ceeds  in  two  stages:  From 


/ 
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dV  (T  (oE^  +  E^^f  >  0 


we  find 


jkj  if^"'"  *  (^''^lih  *  t^''^'  -  ■'w  i° 


Since  this  holds  for  all  a,  it  follows  that 


n    ^\  (2k+l)     (2k+2)         ^         . 

Next,  we  note  that 

JdSn_«E.AH.«  J  dS  E  •   I/.  •     E-   rdSH»^«H<0 

for  any  fields  that  satisfy  the  boimdary  conditioiis.  This  follows  from  the 
interpretation  of  the  Poynting  vector  as  the  flux  of  energy:  for,  since  the 
energy  is  dissipated  only  inside  the  conducting  region,  the  energy  in  the 
external  region  must  have  a  net  flow  in  the  direction  opposite  to  _n  (the  out- 
ward normal  to  S).   Also,  since  the  boundary  conditions  are  derived  by  con- 
sidering the  exact  fields  in  the  external  region,  the  boundary  conditions  will 
insure  that  the  energy  flow  is  in  the  correct  direction. 

We  may  therefore  write,  since  ^^  and  ^^  differ  only  by  a  positive 
constant, 

j  dV  Jv  X  (aEj^  .  E^^^)  r  -  JdS(aEj^  .  ^^^)  .^'  .  (c^  .  E^^^^)  >  0  , 
or,  using  (U.3), 

JdV  <r(aE^  .  E^^^)  .  (aEj^^^  *  5^^)  >  0  . 
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From  this  it  follows  that: 


/.    n\  (2k)  (2k-t-l)    .   „^.  ^  ^ 

(^•7^  T2krn    ^      T2kr-  '   °^-  ^2k    ^      ^2k-.l  • 

In  summary,  we  have 

(h.8)  r^.i  T^,,,  <   T2,,2 <  -r. 

The  remainder  of  the  convergence  proof  proceeds  exactly  as  in  the 
scalar  case:  if  T    is  the  largest  eigenvalue,  3nd  E    is  the  corresponding 
eigenf unction,  then 

^(1)   ^   yd) 

as  k  -»•  CO,  providing  that  the  initial  (generating)  function  is  not  orthogonal 

to  ^   •  A  corresponding  statement  holds  for  the  smaller  eigenvalues  'X       (n  /  1) 

and  their  eigenf iinctions  E.   ,  except  that  now  the  generating  function  is  chosen 

(k) 
orthogonal  to  all  E;   for  k  <  n. 

An  essential  part  of  the  convergence  proof  is  a  statement  concerning 

the  boundedness  of  each  member  of  the  sequence  of  iterated  fTinctions.  For  this 

purpose,  it  is  convenient  to  replace  Eqs,  (U.2)  and  (U.3)  by  an  integral  relation. 

This  may  be  done  by  defining  a  {scalar)   Green's  function  to  satisfy: 

(U.9)  V^   G(r,  r')  -  -5(r  -  r')^ 

(U.IO)  (VG(r,  r')  ;«  I)  X  n  »  Jt*   G(r,  r')        (r  on  S), 

where  I  is  the  unit  dyadic.  We  then  have: 


-  lu- 


(U.ll)  4(r)  -  j  dv'G(r,  r')  ^  <r  (r')  £j^,^(r' ) 

-^  c 

as  the  solution  of  (U.2)  subject  to  the  boundary  condition  (U.3). 


?•  Iteration  of  ^ 

So  far,  we  have  considered  only  the  iteration  of  E,     The  iteration  of 
jj  appears  to  be  more  difficult.  One  method,  however,  is  to  make  the  iteration 
of  K  dependent  on  that  of  E.  by  defining: 

(5.1)  Ay.    -  '  ^  Bj^  , 

where  E,  is  fixed  by  Eqs.  (U.l)-(U.3).  We  now  have  from  (U.2), 


(5.2)  ^^i  T^^y  '  I 


■^-1  • 


The  boundary  condition  on  H,  is  automatically  fixed  by  that  on  Kj  from  (U.3) 
and  (5.1) 


Ji. '»  -^'  -h-y  '  [~f\      '°"''- 


Using  J/  -  cT^'  and  (2.8),  and  defining  V  »  c  ^^,  the  equation 


above  is 
equivalent  to: 


In  Eqs.  (5*2)  and  (U.3),  the  iteration  of  H_ has  been  freed  of  any 
explicit  dependence  on  the  iteration  of  E,  However,  it  appears  from  the  deriva- 
tion of  these  equations  that  the  two  iteration  schemes  are  not  independent,  but 
each  is  equivalent  to  the  other. 
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6,  Mapping 

Two-dimensional  decay  problems  are  amenable  to  treatment  by  conformal 
mapping  methods.  For  example,  taking  Eqs.  (3 •16)  in  two  dimensions,  suppose 
that  we  map  the  z-piane  (z  ■  x  +  iy)  into  the  z  -plane  (z  =  x  +  iy  ): 


(6.1) 


Z  ■  z(z  ) r 


The  line  elements  are  related  by  (we  use  *  for  complex  conjugate); 


2     2       * 
dx  +  dy  ■  dz  dz 


dz 
dz 


(dx  ^  +  dy  ^) 


Therefore,  in  the  z  -plane,  we  get  an  equation  identical  in  form  to  (3.16) 
except  that  o'(z)  is  replaced  by 


(6.2) 


<r    « 


dz 


dz 


<r. 


By  means  of  (6.1)  and  (6.2),  the  conductivities  of  two  problems  having 
the  same  eigenvalue,  X,  may  be  related  to  one  another. 
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